use "C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Bases\Agregadas_Stata\Agregada.dta", clear

estimates drop _all


********************************Modelo para conclusão do EM

*Sem controles
logistic rede_ens b(first).ano##b(1).quintis_renda if (rec_mod==1) & ((ano==1992 | ano==2002 | ano==2022) | (ano==2012 & versao==1)) & (ens_sup_a==1) [pweight = peso] 

*Com controles

*2001
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2001) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2001

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2001

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2001

*2002
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2002) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2002

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2002

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2002

*2003
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2003) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2003

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2003

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2003

*2004
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2004) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2004

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2004

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2004

*2005
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2005) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2005

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2005

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2005

*2006
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2006) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2006

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2006

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2006

*2007
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2007) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2007

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2007

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2007

*2008
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2008) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2008

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2008

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2008

*2009
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2009) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2009

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2009

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2009

*2011
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2011) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2011

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2011

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2011

*2012
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2012) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2012

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2012

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2012

*2013
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2013) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2013

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2013

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2013

*2014
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2014) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2014

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2014

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2014

*2015
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2015) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2015

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2015

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2015

*2016
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2016) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2016

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2016

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2016

*2017
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2017) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2017

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2017

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2017

*2018
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2018) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2018

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2018

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2018

*2019
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2019) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2019

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2019

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2019

*2020
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2020) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2020

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2020

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2020

*2021
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2021) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2021

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2021

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2021

*2022
quietly logistic rede_ens b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2022) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) & (ens_sup_a==1) [pweight = peso] 
estimates store mod_rede_ens_renda_2022

quietly margins quintis_renda, post
estimates store margins_rede_ens_renda_2022

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_rede_ens_renda_2022


****************Tabelas

*OR
etable, estimates(mod_rede_ens_renda_2001 mod_rede_ens_renda_2002 mod_rede_ens_renda_2003 mod_rede_ens_renda_2004 mod_rede_ens_renda_2005 mod_rede_ens_renda_2006 mod_rede_ens_renda_2007 mod_rede_ens_renda_2008 mod_rede_ens_renda_2009 mod_rede_ens_renda_2011 mod_rede_ens_renda_2012 mod_rede_ens_renda_2013 mod_rede_ens_renda_2014 mod_rede_ens_renda_2015 mod_rede_ens_renda_2016 mod_rede_ens_renda_2017 mod_rede_ens_renda_2018 mod_rede_ens_renda_2019 mod_rede_ens_renda_2020 mod_rede_ens_renda_2021 mod_rede_ens_renda_2022) cstat(_r_b, nformat(%4.3f)) stars( , attach(_r_b) ) mstat(N, nformat(%9.0g))

collect layout (coleq#colname#result[_r_b] result[N] result[aic] result[bic] result[chi2] result[df_m] result[p] result[ll] result[r2_p]) (etable_depvar#stars) (), name(ETable)
collect export "C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Resultados\Modelos.xlsx", as(xlsx) sheet(rede_ens_renda_OR, replace) cell(A1) modify


*RR
etable, estimates(nlcom_rede_ens_renda_2001 nlcom_rede_ens_renda_2002 nlcom_rede_ens_renda_2003 nlcom_rede_ens_renda_2004 nlcom_rede_ens_renda_2005 nlcom_rede_ens_renda_2006 nlcom_rede_ens_renda_2007 nlcom_rede_ens_renda_2008 nlcom_rede_ens_renda_2009 nlcom_rede_ens_renda_2011 nlcom_rede_ens_renda_2012 nlcom_rede_ens_renda_2013 nlcom_rede_ens_renda_2014 nlcom_rede_ens_renda_2015 nlcom_rede_ens_renda_2016 nlcom_rede_ens_renda_2017 nlcom_rede_ens_renda_2018 nlcom_rede_ens_renda_2019 nlcom_rede_ens_renda_2020 nlcom_rede_ens_renda_2021 nlcom_rede_ens_renda_2022) cstat(_r_b, nformat(%4.3f)) cstat(_r_lb, nformat(%4.3f)) cstat(_r_ub, nformat(%4.3f)) export("C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Resultados\Modelos.xlsx", as(xlsx) sheet(rede_ens_renda_RR, replace) cell(A1) modify)
